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ABSTRACT 

Natural  convection  liquid  cooling  of  simulated  electronic  compo- 
nents in  a  vertical  channel  was  investigated.  The  test  surface 
contained  a  single  column  of  eight  rectangular,  protruding  heated  ele- 
ments, each  simulating  a  20-pin  dual-in-line  package.  Temperature 
measurements  and  flow  visualization  were  performed  for  a  number  of 
power  dissipation  levels  and  channel  widths.  Collectively,  this  infor- 
mation was  used  in  interpreting  the  flow  and  transport  characteristics. 
A  correlation  to  predict  the  heat  transfer  rates  was  developed  based 
on  the  component  surface  temperatures.  Optimum  channel  widths 
were  determined  from  these  surface  temperature  measurements  for 
the  range  of  power  levels  investigated.  Temperature  distributions  in 
the  fluid  were  measured  using  a  traversing  thermocouple  probe. 
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I.  INTRODUCTION 

A.   STATEMENT  OF  PROBLEM 

The  continuing  trend  toward  microminiaturization  of  electronic 
components  has  made  thermal  design  critical  to  achieving  higher 
packaging  densities.  The  ever- increasing  heat  fluxes  must  be  removed 
while  keeping  the  device  junction  temperatures  typically  below  100°C 
[Ref.  1].  The  failure  of  microelectronic  chips  increases  exponentially 
as  the  junction  temperatures  increase  [Ref.  2]. 

Microchips  act  as  discrete  sources  of  thermal  energy.  Effective 
removal  of  the  generated  energy  is  of  major  import  in  determining 
how  they  might  be  used  without  overheating  the  elements  themselves 
or  the  surrounding  package.  Thermal  control  is  one  of  the  critical 
considerations  in  the  design  of  any  electronic  cooling  system.  The 
thermal  designer's  task  is  further  complicated  by  the  presence  of 
arrays  of  electronic  components  placed  in  proximity  to  each  other, 
requiring  detailed  considerations  of  transport  interactions  between 
neighboring  elements. 

The  design  of  packaged  electronic  components  must  consider  all 
anticipated  modes  of  heat  transfer.  A  complete  thermal  analysis  will 
account  for  all  three  mechanisms  of  thermal  energy  transport: 

•  conduction 

•  convection 

•  radiation 


The  areas  of  concern  include  transferring  heat,  via  conduction,  from 
internal  heat  sources  to  the  external  surface  of  the  package  and  the 
subsequent  removal  of  heat  from  the  external  surface  by  a  cooling 
medium.  Also,  the  coolant  temperature  must  be  maintained  at  pre- 
scribed levels  for  a  given  environment. 

The  present  work  focuses  on  the  heat  removal  from  electronic 
packages  through  convective  cooling.  This  convection  may,  in  general, 
be  forced  or  natural;  the  medium  liquid  or  gas.  Because  of  the  higher 
density  and  thermal  conductivity  of  liquids  compared  to  air,  heat 
transfer  coefficients  an  order  of  magnitude  or  more  larger  than  for  air 
are  attainable.  Due  to  the  significantly  high  heat  transfer  rates  achiev- 
able, cooling  of  electronic  components  by  direct  liquid  immersion  has 
been  identified  as  an  extremely  effective  thermal  control  technique 
[Ref.  2]. 

Convection  heat  transfer  rate  from  a  surface  at  temperature  Ts  of 
area  A  to  surrounding  fluid  environment  Tj^f  is  given  by  Newton's  law 
of  cooling: 

q  =  Ah{Ts  -  Tm) 

This  is  a  definition  of  h,  the  convective  heat  transfer  coefficient,  more 
than  a  phenomenological  law.  Evaluation  of  h  depends  upon  heat 
source  geometry,  orientation,  fluid  properties,  and  fluid  velocity. 

Natural  convection  results  when  fluid  density  gradients  exist  in 
the  presence  of  a  body  force  field  such  as  gravity.  Density  gradients  in 
the  fluid  commonly  occur  due  to  temperature  gradients.   Advantages  of 


employing  natural  convection  in  electronic  cooling  include  high  relia- 
bility, low  cost,  reduced  noise,  and  a  simpler  design.  Despite  these 
potential  benefits,  only  a  limited  number  of  studies  are  available  on 
liquid  immersion  natural  convection  cooling  of  electronic  components. 

Most  existing  natural  convection  studies,  both  experimental  and 
analytical,  have  considered  heating  surfaces  that  were  semi-infinite. 
Electronic  component  cooling  configurations  typically  employ  arrays  of 
small,  discrete,  flush,  or  protruding  heat  sources.  The  available  heat 
transfer  correlations  are  inadequate  in  predicting  transport  in  these 
geometries.  Convective  heat  transfer  from  an  array  of  microchips  is 
strongly  dependent  on  package  shape  and  array  density.  Indeed,  in 
real  systems,  individual  elements  or  modules  often  vary  significantly  in 
size,  shape,  and  in  the  manner  in  which  they  are  arranged  on  the 
substrate. 

The  geometric  complexity  is  further  compounded  by  the  nonuni- 
form nature  of  heat  dissipation  from  the  array.  Individual  component 
dissipation  is  a  function  of  its  utilization.  Variation  in  power  dissipa- 
tion from  component  to  component  on  the  same  board  and  from  board 
to  board  within  the  system  usually  exists.  The  orientation  of  the 
boards  on  which  the  elements  are  mounted  is  yet  another  variable. 

Often,  these  semi-regular  arrays  face  the  back  of  an  adjacent 
board.  The  resulting  vertical  channel  has  one  surface  which  is  rela- 
tively smooth  and  the  other  one  covered  by  an  arrangement  of  large, 
heat-dissipating  elements.  Within  the  channel,  local  convective  heat 
transfer  is  driven  by  two  phenomena,  a  local  buoyancy  force  due  to  the 


heated  components  and  a  forced  convection  effect  arising  from  the 
flow  history,  over  the  entire  channel  length,  including  both  heated  and 
unheated  regions.  A  simultaneous  consideration  of  these  mechanisms 
is  essential  in  obtaining  predictive  capability  for  transport  rates  in 
these  systems. 

The  complexity  of  the  problem  described  demands  an  approach 
which  builds  on  the  fundamental  physical  processes.  Efficient  cooling 
cannot  be  attained  without  understanding  the  heat  transfer  from  each 
package  by  a  determination  of  its  flow  and  thermal  fields.  The  overall 
goal  is  to  develop  a  methodology  based  on  sound  physical  principles 
which  will  allow  the  prediction  of  individual  component  temperatures, 
by  superposition,  in  a  nonuniformly  heated  array  of  components. 

B.      IMMERSION    COOLING— ANALYTICAL    AND    EXPERIMENTAL 
STUDIES 

Reference  3  provides  an  excellent  literature  review  of  previous 
work,  both  analytical  and  experimental,  in  the  field  of  liquid  immer- 
sion cooling  of  electronic  equipment.  Additionally,  very  relevant  work 
has  been  performed  by  Ortega  and  Moffat  [Refs.  4-6].  They  have 
experimentally  examined  natural  convection  air  cooling  of  simulated 
electronic  components  protruding  from  a  vertical  surface,  with  and 
without  a  shrouding  wall.  Heat  transfer  rates  significantly  higher  than 
for  a  smooth  vertical  plate  were  found.  They  interpreted  this  as  an 
indication  of  turbulent  free  convection.  When  a  shroud  was  in  place, 
they  found  plate -averaged  heat  transfer  coefficients  40  to  50  percent 
higher  than  for  a  smooth,  parallel  plane  channel. 


Jaluria  [Ref.  7]  performed  two-dimensional  computations  for 
natural  convection  air  flow  over  a  vertical  surface  with  multiple  wide, 
flush-mounted  heaters.  The  study  treated  the  flow  as  a  boundary  layer 
problem.  Rajakumar  and  Johnson  [Ref.  8]  reported  computations  of 
free  convection  heat  transfer  using  rectangular  strip  heating  surfaces. 
Their  finite  element  approach  showed  significant  stagnant  areas 
between  the  blocks. 

C.     OBJECTIVE 

An  experimental  study  was  undertaken  to  to  determine  the 
hydrodynamic  and  thermal  characteristics  of  buoyancy-induced  flow 
over  several  heated  protruding  components  mounted  on  a  vertical 
surface. 

A  single  column  array  was  chosen  because  it  provided  the  sim- 
plest geometry  featuring  discretized  heat  dissipation.  The  effects  of  a 
neighboring  card  were  simulated  by  a  movable  vertical  shroud  plate. 
Specific  objectives  of  this  study  were: 

•  to  visualize  the  transient  and  steady-state,  natural  convection  flow 
within  an  interrupted  channel,  for  a  range  of  component  power 
dissipation  levels  and  channel  widths; 

•  to  measure  component  temperatures  for  various  power  inputs  and 
shroud  spacings  to  develop  an  appropriate  nondimensional  cor- 
relation to  predict  the  heat  transfer  rates  for  this  geometry;  and 

•  to  measure  temperature  distributions  in  the  adjacent  fluid,  using  a 
traversing  thermocouple  probe  to  assist  in  evaluating  the  hydro- 
dynamic  and  thermal  nature  of  the  flow. 


This  study  was  a  follow-on  to  work  described  in  Reference  3  and 
is  part  of  a  continuing  effort  to  study  natural  convection  liquid  immer- 
sion cooling  of  electronic  equipment. 


II.    EXPERIMENT 

A.     GENERAL  DESIGN  CONSIDERATIONS 

Detailed  descriptions  of  equipment  design  and  construction  are 
available  in  Reference  3.    Only  a  brief  overview  is  presented  here. 

As  seen  in  Figure  1 ,  the  assembly  consisted  of  a  vertical  test  sur- 
face which  contained  a  single  column  of  eight  rectangular  protruding 
heated  elements.  Each  protrusion  simulated  a  20-pin  dual-in-line 
package.  The  elements  were  stainless  steel  blocks  spaced  one  inch 
apart,  center  to  center,  on  a  plexiglass  substrate.  A  12.70  mm  thick 
plexiglass  shroud,  when  used,  simulated  the  reverse  side  of  a  second 
printed  circuit  board. 

The  blocks  (Figures  2  and  3)  had  imbedded  thermocouples  on 
each  of  the  five  surfaces  in  contact  with  the  fluid  to  provide  surface 
temperature  measurements.  A  sixth  thermocouple  was  placed  in  the 
center  of  each  element  mounting  slot  for  measurements  of  heater 
temperature.  A  miniature  thermocouple  probe  mounted  on  a  three- 
dimensional  traverse  enabled  measurements  of  temperature  distribu- 
tions within  the  fluid. 

A  foil  heater  mounted  on  the  component  side  facing  the  substrate 
was  powered  using  a  regulated  D.C.  power  supply.  Each  heater  was 
run  in  series  with  a  precision  resistor  and  all  eight  heaters  were  in 
parallel  with  the  power  supply.     The  current  to  each  heater  was 
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Figure  1.    System  Configuration 
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Figure  2.   Heater  Block  Schematic 

(after  Hazard,  Ref.  3) 
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Figure  3.    Mounted  Foil  Heater.  Two  Perspectives 

(after  Hazard,  Ref.  3) 
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calculated  by  subtracting  the  heater  voltage  from  the  source  voltage, 
then  dividing  by  the  resistance  of  the  precision  resistor.  Power  dis- 
sipation through  each  element  was  determined  by  calculating  the 
product  of  heater  voltage  and  heater  current. 

Figure  2  shows  the  orientation  of  the  coordinate  axes  used  in  the 
discussion  to  follow.  The  front  face  of  the  component  is  that  surface 
in  the  x-z  plane  in  contact  with  the  fluid.  The  side  faces  are  the  two 
x-y  surfaces,  right  and  left  as  determined  by  viewing  the  component 
from  the  front.  The  top  and  bottom  faces  are  those  in  the  y-z  plane, 
downstream  and  upstream  of  buoyant  fluid  flow,  respectively. 
Component  height  is  in  the  x-direction,  depth  in  y,  and  breadth  in  the 
z-direction. 

The  immersion  fluid  was  purified  water  contained  in  a  one  cubic 
meter,  plate  glass  tank.  Suspended  in  the  water  were  particles  of 
Pliolite,  an  inert  pigment  with  a  specific  gravity  of  0.93.  The  particles 
were  illuminated  with  an  eight  milliwatt  helium  neon  laser  split  into  a 
plane  by  use  of  a  cylindrical  lens.  The  plane  of  visualization  was  varied 
by  realignment  of  the  laser  lens  assembly  and  the  camera. 

B.     EXPERIMENTAL  PROCEDURE 

1.  Established  Method 

Details  are  again  provided  in  Reference  3,  hereafter  referred 
to  as  the  previous  study.    Noteworthy  differences  are  noted  here. 

2.  Instrument  Settings 

Power  settings  of  0.2,  0.5,  1.0,  1.2,  and  1.5  watts  were  used 
for  each  channel  width.     Because  of  the  inability  to  de-aerate  the 
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immersion  fluid,    1.5  watts  was  the  upper  limit  of  power  available. 
Settings  above  this  generated  bubbles  on  the  block  surfaces. 

3.  Instrument  Readings 

When  measuring  steady-state  surface  temperatures,  succes- 
sive temperature  measurements  at  10-minute  intervals  were  com- 
pared. Only  when  all  thermocouples  varied  less  than  0.10°  C  were  the 
final  measurements  taken.  The  temperature  acquisition  program  of 
the  previous  study,  modified  to  reflect  the  numbering  of  elements 
from  bottom  to  top  (flow  direction)  and  to  allow  data  storage 
(Appendix  A),  was  then  executed. 

When  using  the  traversing  thermocouple  probe,  the  data 
acquisition  system  was  programmed  to  continuously  monitor  the  volt- 
age registered  by  a  single  channel.  Observation  of  the  on-screen  print- 
out allowed  determination  of  temperature  variation  from  point  to  point 
within  the  flow  field,  as  well  as  at  any  point  over  time. 

A  program  (Appendix  A)  which  monitored  from  one  to  eight 
thermocouples  over  time  was  used  when  making  a  record  of  the  tran- 
sient temperature  response  of  the  surface  thermocouples  to  power 
step  up  or  step  down.  Typical  data  presented  were  for  the  front  face 
temperatures  of  the  various  components. 

4.  Photographic  Technique 

Best  results  were  obtained  when  using  an  f  stop  of  2.8  and 
exposure  time  of  20  to  40  seconds.  Camera  placement  was  approxi- 
mately six  inches  from  the  tank  wall  at  a  height  which  allowed  all 
eight  blocks  in  the  field  of  view.    The  camera  was  inclined  at  an  angle 
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of  about  5  degrees  from  the  test  surface  when  a  shroud  was  present. 
This  was  helpful  in  minimizing  reflected  light.  A  shroud  spacing  of 
two  times  component  height  from  the  plexiglass  substrate  was  deter- 
mined to  be  the  minimum  at  which  useful  pictures  could  be  obtained. 

C.     DATA  ANALYSIS 

While  Qcond  and  Qconv  were  calculated  exactly  as  recommended 
by  the  previous  study,  further  data  reduction  differed  significantly. 

For  each  block,  a  surface -averaged  temperature  was  generated  by 
multiplying  the  temperature  and  area  of  each  fluid -exposed  face  and 
dividing  the  summation  over  the  block  by  the  total  exposed  face  area 
of  the  block.  The  average  over  the  eight  blocks  of  the  surface-averaged 
temperatures  was  then  averaged  with  the  immersion  fluid  tempera- 
ture to  produce  a  film  temperature.  Fluid  properties  based  on  curve 
fits  (Appendix  B)  of  data  from  Incropera  and  Dewitt  [Ref.  9]  were  then 
calculated  for  the  film  temperature. 

A  new  characteristic  length,    L  ,   for  determining  a  nondimen- 

Q 

sional  temperature  excess,  77-,  a  modified  Grashof  number,  Gr*,  and 
Nusselt  number,  Nu,  was  chosen.  The  progression  of  its  use  for  each 
block  is  as  follows: 
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The  various  symbols  used  are  defined  in  the  List  of  Symbols  section. 

Other  characteristic  lengths,  L  to  determine  block  spacing  along 
the  plexiglass  substrate  and  D  to  represent  the  horizontal  breadth  of 
the  element,  were  also  chosen.    Their  use  is  as  follows: 
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-y-    =    nondimensional  location  of  component  centers  along  the 
test  surface 


-J-    -    horizontal  position  along  the  element's  front  face 

Quantitative  development  was  best  expressed  as  plots  of: 

•  Nuvs.  Gr*  relationships 

•  nondimensional  block  surface  temperature  excess  vs.   element 
position 

•  temperature  distributions  in  the  fluid  adjacent  to  the  test  surface 

•  block  surface  temperatures  vs.  time 

Pertinent  trends  in  these  variations  are  discussed  next. 
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HI.    NO  SHROUD- INFINITE  CHANNEL  WIDTH 

A.     ALL  ELEMENTS  POWERED,  STEADY  STATE 
1.     Flow  Visualization 

Photographs  of  the  natural  convection  flow  in  a  plane  (x-y) 
geometrically  centered  through  each  component  are  shown  in  Fig- 
ure 4.  The  visualizations  represent  all  elements  heated  to  power 
settings  of  0.2,  0.5,  and  1.0  watts,  left  to  right,  respectively.  The 
exposure  times  were  40,  20,  and  20  seconds,  left  to  right.  The  pho- 
tograph toning  techniques  required  for  mass  printing  of  this  work  do 
not  allow  sufficient  resolution  for  detailed  study.  Should  the  reader 
desire  to  see  prints  of  much  greater  clarity,  contact  Y.  Joshi  at  the 
address  given  in  the  distribution  list  of  this  study. 

The  flow  appears  to  consist,  in  each  case,  of  strong  upward 
motion  driven  by  the  buoyant  force  near  the  test  surface.  The  fluid 
tends  to  follow  the  contours  of  the  protrusions.  The  particle  traces 
are  closer  passing  the  element  front  face  as  compared  to  the  inter- 
block space.  Dead  zones  are  not  obvious  except  perhaps  for  very  small 
areas  where  the  blocks'  top  faces  intersects  the  substrate  and  for  a 
slightly  larger  area  above  the  top  face  of  the  last  block  in  the  array. 
There  is  no  evidence  of  recirculation  or  formation  of  dead  zones  near 
the  bottom  face  of  any  block. 
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Figure  4.    Steady  Flow  in  X-Y  Plane  at  Power  Levels  of  0.2,  0.5, 
and  1.0  Watts.    Powers  Increase  From  Left  to  Right 
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In  all  three  photographs,  a  developing  boundary  layer  characteristic  is 
evident  farther  out  from  the  protrusions.  It  attains  a  thickness  of 
about  3.5  component  depths  for  0.2  watts  and  2  for  1.0  watt.  Its 
origin  moves  upstream  with  increased  power,  from  about  1  com- 
ponent height  at  0.2  watts  to  3  component  heights  at  1.0  watt.  The 
velocities  appear  greater  and  the  boundary  layer  becomes  thinner  as 
power  is  increased.  Ar  3nt  fluid  is  entrained  over  the  entire  length 
of  the  layer. 

Figure  5  presents  flow  visualization  in  the  plane  (x-z)  passing 
through  the  front  faces  of  the  elements.  All  elements  were  powered 
to  1.2  watts  and  exposure  time  was  40  seconds.  The  photograph 
shows  an  upward-moving  buoyant  layer  near  the  components.  Ambient 
fluid  is  entrained  from  each  side  into  the  flow,  resulting  in  an  increase 
in  the  horizontal  extent  of  the  upflow  region  dowr,  stream.  The  flow  in 
the  vicinity  of  the  elements  appears  three-dimensional. 

As  with  the  sid  view,  the  flow  appears  to  follow  the  contours 
of  the  projections  to  about  the  same  degree  (dipping  in  about  3-5  mm 
from  the  side  faces),  the  flow  also  appears  to  wrap  around  the  cor- 
ners. The  particle  traces  about  20  percent  along  the  block  (z  direc- 
tion) from  the  side  edge  appear  to  be  largely  in  the  x-y  plane  and 
those  at  the  side  face  are  in  the  x-z  plane.  The  particle  traces  in 
between  are  3-D,  each  progressively  moving  along  the  90°  arc. 

Overall,  the  flow  can  best  be  described  as  following  the  con- 
tour of  an  hourglass  split  in  the  x-z  plane.   The  elements  are  the  upper 
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Figure  5.    Flow  in  X-Z  Plane.  1.2  Watts 
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and  lower  vessels.    This  configuration  is  repeated  over  the  length  of 
the  array.     It  becomes  more  pronounced  as  the  flow  moves  up  the 
array  and  is  eventually  shed  as  a  plume  above  the  last  element. 
2.      Quantitative 

Figures  6  and  7  are  nondimensional  representations  of  tem- 
perature excess  for  0.2  and  1.0  watts,  measured  on  various  component 
surfaces  vs.  block  position.  Also  shown  are  the  data  for  heater 
temperature  measurements.  The  two  plots  follow  the  same  patterns. 
The  numbers  for  the  0.2  watts  are,  as  expected,  larger  due  to  the 
inverse  relationship  with  T*,  which  is  h  dependent.  Also  as  antici- 
pated, the  heater  excesses  are  much  greater.  The  drop  in  the  heater 
temperature  between  blocks  6  and  7  could  perhaps  result  from  a 
slightly  smaller  (-4%  less)  power  achieved  by  block  7.  Unfortunately, 
the  heater  thermocouple  for  block  8  was  inoperative.  For  both  power 
settings,  the  imbedded  surface  thermocouples  showed  increasingly 
greater  excesses  along  the  array. 

Of  primary  importance  in  thermal  design  is  the  availability  of 
the  appropriate  heat  transfer  relationships.  As  discussed  earlier,  such 
correlations  for  discrete  protruding  heat  sources  in  liquids  have  not 
been  available  in  literature  up  to  now.  To  this  end,  a  logarithmic  plot 
of  Nu  vs.  Gr*  (Figure  8)  was  generated  with  eight  curves  representing 
element  position.  Each  curve  had  five  points  for  the  power  settings 
0.2,  0.5,  1.0,  1.2,  and  1.5  watts  resulting  in  a  range  of  Gr*.  Although  a 
linear  pattern  with  Gr*  was  evident,  the  spread  in  data  indicated  an 
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Figure  6.   Nondimensional  Temperature  Excess  vs. 
Position.  0.2  Watts 
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Figure  7.    Nondimensional  Temperature  Excess  vs. 
Position.  1.0  Watt 
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Figure  8.    Data  for  Nusselt  Number  vs. 
Flux-Based  Grashof  Number 
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additional  dependence  on  element  location.  A  review  of  the  data  sug- 
gested a  method  by  which  it  could  be  collapsed.  A  modified  Nusselt 
number,  Nu*,  was  calculated  by  use  of  a  scaling  factor  based  on  the  net 
convected  energy  of  the  flow  at  a  downstream  location.  Since  each 
block  dissipates  approximately  an  equal  amount  of  power,  the  net 
convected  energy  is  proportional  to  the  number  of  upstream  compo- 
nents. The  scaling  factor  was  of  the  form  Bna  where  Bn  is  the  block 
number  and  a  >  0.  A  value  of  a  =  1/6  was  found  to  provide  the  best 
correlation  for  the  data.  As  expected,  the  strong  three-dimensional 
effect  on  transport  results  in  a  weaker  downstream  dependence  than 
the  1/5  exponent  in  the  Fujii  and  Fujii  correlation  for  a  uniform  flux, 
semi-infinite  vertical  surface.    Therefore, 

Nu*=  Nu  ■  Bn6 

for  each  element  at  each  power  setting.  The  modified  data  are  plotted 
in  Figure  9  along  with  the  linear  least  squares  fit.  Based  on  this,  the 
following  correlating  equation  for  the  heat  transfer  rates  is  obtained: 


Nu*=  L88(Gr*)15 


for  50  <  Gr*  <  1000  and  Pr  ~  7. 
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Figure  9.    Nusselt  Number  vs.  Flux-Based  Grashof  Number 

With  Least  Squares  Fit 
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B.     ALL  ELEMENTS  POWERED,  TRANSIENT 

1.  Flow  Visualization 

Figure  10  is  a  series  of  photographs  taken  during  the  0.2  watt 
step-up  transient.  For  each  picture,  the  exposure  time  was  20  sec- 
onds. Exposure  of  the  first  picture  began  with  power  application.  The 
others  were  initiated  at  20,  80,  120,  and  200  seconds,  respectively, 
following  the  start  of  heating.  Clearly,  the  first  two  show  predom- 
inantly parallel  flow.  The  next  pictures  represent  increasing  devel- 
opment of  the  contour  flow  and  boundary  layer  characteristic 
development.  By  200  seconds,  almost  fully  developed  flow  is 
recognized. 

Figure  1 1  is  a  similar  series,  again  of  20-second  exposures, 
taken  during  the  step-down  transient  from  1.0  watt.  The  first  visual- 
ization at  far  left  began  at  power  down.  The  remaining  photographs 
were  initiated  80,  160,  240,  and  440  seconds,  respectively,  following 
shutdown.  The  steady-state  flow  characteristics  persist  to  an 
increasingly  lesser  degree  for  the  entire  period  observed.  The 
entrainment  of  ambient  fluid  lessens  at  large  times.  Additionally,  flow 
origin  moves  slightly  downstream  over  time. 

2.  Quantitative 

The  1.0-watt  step-up,  step-down  transient  temperature 
responses  for  the  component  front  surfaces  are  presented  on  the 
same  graph  (Figure  12).  Steady  state  is  achieved  about  25  percent 
faster  for  the  power  up.    The  curves  representing  blocks  1  through  8 
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Figure  10.    0.2  Watt  Step-Up  Transient  During 
Various  Time  Intervals 
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Figure  11.    1.0  Watt  Step-Down  Transient  During 
Various  Time  Intervals 
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Figure  12.   Transient  Temperature  Response  Following 
Initiation  of  Heating.  1.0  Watt 


29 


excluding  3  travel  closely  together,  especially  in  the  first  quarter  of 
the  response  period  for  the  step.  The  almost  single  line  corresponds 
to  the  parallel  flow  pattern  seen  in  the  photographs.  This  indicates  a 
purely  diffusive  transport  response  for  a  large  part  of  the  transient. 
Later,  convective  effects  become  appreciable  as  entrainment  develops 
downstream.  For  the  step  down,  steady  state  temperature  is  reached 
much  more  quickly  than  breakdown  of  the  established  flow.  The  per- 
sistence of  flow  is  conducive  to  the  rapid  cooling  of  the  blocks. 

The  plot  for  0.2  watts  (Figure  13)  reaches  steady  state  at 
about  the  same  time  for  both  step  up  and  step  down.  This  period  is 
longer  than  for  either  1.0-watt  case  and  indicates  that  the  stronger 
flow  more  readily  establishes  the  final  thermal  pattern. 

C.     TEMPERATURE  MEASUREMENTS  IN  THE  FLUID 

Characteristics  of  the  thermal  transport  in  the  fluid  were  mea- 
sured using  a  traversing  thermocouple  probe  adjacent  to  all  protruding 
elements  and  for  points  midway  between  the  components  (Figures  14 
and  15).  The  fluid  temperature  responses  near  the  blocks  and  inter- 
block spaces  follow  the  same  pattern,  with  temperature  levels  in  the 
vicinity  of  blocks  4  and  6  surprisingly  greater  than  for  the  others. 
Indeed,  both  graphs  can  be  seen  to  have  four  groupings.  Block  1  by 
itself,  the  lowest  curve,  is  followed  by  2,  3,  and  5;  then  4,  7,  and  8; 
and  finally  block  6  by  itself.  More  detailed  local  measurements  are 
needed  to  explain  these  variations. 
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Figure  13.   Transient  Temperature  Response  Following 
Termination  of  Power.  0.2  Watts 
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Figure  14.    Fluid  Temperature  Distribution  in  the  Normal 
Direction.  Awav  From  Block  Faces 
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Figure  15.    Fluid  Temperature  Distribution  in  the  Normal 
Direction.  Awav  From  Plexiglass  Substrate 


33 


The  layer  between  the  blocks  is  slightly  thicker  with  respect  to 
that  at  the  block  faces.  This  might  be  expected  from  the  flow 
visualization  (Figure  4)  showing  that  the  particle  traces  spread  out 
between  elements.  In  each  case,  the  thermal  layer  is  about  one 
component  height  or  less  and  is  therefore  deeply  imbedded  within 
the  momentum  exchange  region.  This  is  expected  due  to  the 
relatively  large  Pr  ~  7  for  water. 

Next,  the  traversing  probe  was  placed  in  the  thermal  layer,  1  mm 
from  the  surface,  and  moved  in  the  x-z  plane  across  the  test  surface 
face  (Figures  16  and  17).  The  temperature  variations,  with  position, 
reinforce  the  photographs  in  suggesting  strongly  three-dimensional 
transport.  Ambient  temperature  is  the  lowest  point  shown.  The 
peaks  at  about  15  percent  from  each  edge  (Figure  16)  correlate  with 
the  regions  in  the  visualization  (Figure  5),  where  the  flow  begins  to 
shift  from  x-y  dominant  to  x-z  and  may  indicate  a  common  region 
through  which  those  particles  are  flowing.  The  result  is  more, 
warmer  fluid  convected  over  that  region.  Likewise,  the  peaks  of  Fig- 
ure 17  might  represent  a  similar  flow  shift  at  locations  between  the 
components. 

Finally,  the  probe  was  placed  above  the  center  of  the  top  face  of 
block  8,  2  mm  out  from  the  substrate  to  measure  the  downstream 
temperature  decay  (Figure  18).  The  extremely  high  temperature  in 
the  immediate  vicinity  of  the  top  face  suggests  a  dead  zone,  as  did  the 
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Figure  16.    Spanwise  Fluid  Temperature  Distribution 
Across  Components 
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Figure  17.    Spanwise  Fluid  Temperature  Distribution 
Across  Intercomponent  Spaces 


36 


Figure  18.    Fluid  Temperature  Decay  Downstream 
of  Top  Component 
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visualizations,  where  a  small  pool  of  fluid  is  heated  primarily  by  ther- 
mal conduction  from  the  component  upper  surface.  There  was  no 
sign  of  a  vortex. 

All  measurements  with  the  probe  were  remarkably  consistent 
over  time  (±1  p.v)  at  any  position.  No  evidence  by  any  measurement 
indicated  the  presence  of  turbulent  flow  for  the  power  levels 
investigated. 

D.     INDIVIDUAL  ELEMENTS  POWERED 

Individual  blocks  were  powered  to  1.0  watt,  one  at  a  time.  The 
resulting  steady-state  array  temperatures  were  monitored  and  the  flow 
visualized.    These  are  discussed  next. 

1.      Flow  Visualization 

A  series  of  photographs  show  single-element  heating  as  that 
element  changes  from  one  block  to  the  next  up  the  array.  The  expo- 
sure times  for  all  pictures  presented  were  40  seconds.  Figure  19  pre- 
sents left  to  right  blocks  1,  3,  and  4  heated  and  Figure  20  shows 
blocks  5,  6,  and  7.  The  flow  with  only  the  lowest  block  heated  resem- 
bles the  all-components-heated  case,  indicating  flow  around  block 
contours  and  a  boundary  layer  characteristic  over  the  entire  length.  As 
progressively  higher  blocks  are  heated,  a  quiescent  area  begins  to 
develop  adjacent  to  the  lower  part  of  the  test  surface.  Typically,  fluid 
entrainment  results  in  vigorous  flow  around  the  first  block  upstream  of 
the  one  heated.  At  two  blocks  upstream,  a  comparatively  weaker  flow 
is  seen  and  at  three  the  fluid  is  nearly  quiescent. 
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Figure  19.   Steady  Flow  Responses  With  Blocks 
1.  3.  and  4  Individually  Heated 
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Figure  20.    Steady  Flow  Responses  With  Blocks 
5.  6.  and  7  Individually  Heated 
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2.      Quantitative 

Nondimensionalized  temperature  excesses  for  the  faces  and 
heater  of  each  element  are  displayed  (Figures  21  through  25).  The 
plots  show  very  similar  temperature  excess  patterns  around  the 
heated  element  as  the  selected  component  progresses  up  the  array. 
The  protruding  elements  have  minimal  effect  when  they  are  part  of 
the  unheated  starting  length. 

E.      COMBINATIONS  OF  ELEMENTS  HEATED 

As  in  Section  D,  progression  was  made  up  the  array,  except  that 
all  upstream  blocks  remained  heated.  Figures  26  through  30  indicate 
very  little  effect  of  powering  of  downstream  components  on  the  tem- 
perature responses  of  already-powered  upstream  components.  Also 
evident  is  a  very  sharp  decay  in  temperature  in  the  unheated  blocks. 
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Figure  21.    Nondimensional  Temperature  Excess  Levels  at 
Front  Faces  of  Individually  Heated  Components 
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Figure  22.    Nondimensional  Temperature  Excess  Levels  at 
Right  Faces  of  Individually  Heated  Components 
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Figure  23.    Nondimensional  Temperature  Excess  Levels  at 
Top  Faces  of  Individually  Heated  Components 
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Figure  24.    Nondimensional  Temperature  Excess  Levels  at 
Bottom  Faces  of  Individually  Heated  Components 
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Figure  25.    Nondimensional  Temperature  Excess  Levels  at 
Heaters  of  Individually  Heated  Components 
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Figure  26.    Nondimensional  Temperature  Excess  Levels  at 
Front  Faces  of  Multiple  Heated  Components 
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Figure  27.    Nondimensional  Temperature  Excess  Levels  at 
Right  Faces  of  Multiple  Heated  Components 
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Figure  28.    Nondimensional  Temperature  Excess  Levels  at 
Tod  Faces  of  Multiple  Heated  Components 
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Figure  29.    Nondimensional  Temperature  Excess  Levels  at 
Bottom  Faces  of  Multiple  Heated  Components 
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Figure  30.   Nondimensional  Temperature  Excess  Levels  at 
Heaters  of  Multiple  Heated  Components 
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IV.    SHROUD  IN  POSITION— FINITE  CHANNEL  WIDTHS 

A  plexiglass  plate  (shroud)  was  used  (Figure  1)  to  simulate  the 
channel  formed  by  two  adjacent  parallel  printed  circuit  boards.  The 
test  surface  as  before  has  protruding  heated  elements.  The  shroud 
simulates  the  reverse  side  of  the  second  board  in  an  array. 

A.     ALL  ELEMENTS  POWERED,  STEADY  STATE 
1.     Flow  Visualization 

Photographs  presenting  natural  convection  flow  in  the  x-y 
plane,  geometrically  centered  through  each  component  (Figure  31), 
were  evaluated.  The  visualizations  are  for  all  elements  heated  to  0.5, 
1.0,  1.2,  and  1.5  watts.  The  exposure  times  were  20  seconds.  The 
shroud  distance  was  18  mm  from  the  substrate  or  approximately  three 
component  depths.  The  first  six  blocks  are  clearly  visible  in  each  pic- 
ture. Figure  32  shows  the  top  portion  of  the  test  surface  for  0.5  watts. 
The  other  power  settings  do  not  differ  significantly  and  are  not  shown. 

The  flow  near  the  elements  appears  similar  to  that  observed 
in  the  absence  of  the  shroud.  However,  the  shrouding  wall  prevents 
ambient  fluid  from  entrainment  into  the  buoyant  up-flow  in  this  plane. 
As  power  is  increased,  the  paths  described  by  the  particle  traces  are 
more  closely  spaced,  indicative  of  increased  velocity  and  increased 
fluid  flow. 
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Figure  31.    Steady  Flow  in  X-Y  Plane.  0.5.  1.0.  1.2.  and  1.5  Watts. 
18  mm  Channel  Spacing.  Bottom  Blocks 
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Figure  32.    Steady  Flow  in  X-Y  Plane.  0.5  Watts.  18  mm 
Channel  Spacing.  Top  Blocks 
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Figure  33  is  the  front  view  of  flow  in  the  x-z  plane  through 
the  front  surface  of  the  elements.  It  is  for  1.0  watt,  18  mm  shroud 
spacing,  and  exposure  time  of  40  seconds.  Unlike  the  no-shroud  con- 
dition, flow  is  more  restricted  to  the  sides  of  the  protruding  blocks. 
The  characteristic  hourglass  flow  seen  for  infinite  channel  width  is 
weak  in  the  vicinity  of  the  middle  two  blocks  and  absent  elsewhere. 
The  entrained  fluid  turns  very  sharply  into  the  flow  direction.  Because 
the  flow  is  restricted  to  the  channel  width  and  is  predominant  on  the 
block  sides,  its  velocity  is  much  greater  here.  Consequently,  there  is 
more  forceful  entrainment  of  fluid  in  the  x-z  plane. 
2.     Quantitative 

Due  to  the  large  amount  of  data  collected,  temperature 
excesses  for  power  settings  of  1.2  and  0.2  watts  only  are  presented  in 
graphical  form.  These  settings  are  representative  of  the  upper  and 
lower  ends  of  the  power  range  available.  Channel  widths  presented 
are  6.2,  8,  14,  18,  and  23  mm  measured  from  the  substrate. 

The  first  width  was  just  barely  greater  than  component 
height.  This  spacing  is  representative  of  pure  conduction  dominated 
transport.  The  nondimensionalized  temperature  excess  vs.  element 
position  graphs  for  0.2  watts  (Figures  34-38)  and  1.2  watts  (Figures 
39-43)  show  significant  difference  between  the  smallest  spacing  and 
all  others.  With  that  spacing  excluded,  Figures  44  through  53  present 
the  relevant  data  in  a  manner  better  suited  for  analysis.  The  plots  for 
0.2  watts  show  the  next  closest  spacing,  8  mm  from  the  substrate,  to 
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Figure  33.    Flow  in  X-Z  Plane.  1.0  Watt.  18  mm  Channel  Spacing 
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Figure  34.    Nondimensional  Temperature  Excess.  0.2  Watts. 
Various  Spacings.  Front  Faces 
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Figure  35.    Nondimensional  Temperature  Excess.  0.2  Watts. 
Various  Spacings.  Right  Faces 
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Figure  36.    Nondimensional  Temperature  Excess.  0.2  Watts. 
Various  Spacings.  Top  Faces 
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Figure  37.    Nondimensional  Temperature  Excess.  0.2  Watts. 
Various  Spacings.  Bottom  Faces 
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Figure  38.    Nondimensional  Temperature  Excess.  0.2  Watts. 

Various  Spacings.  Heaters 
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Figure  39.    Nondimensional  Temperature  Excess.  1.2  Watts. 
Various  Spacings.  Front  Faces 
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Figure  40.    Nondimensional  Temperature  Excess.  1.2  Watts. 
Various  Spacings.  Right  Faces 
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Figure  41.    Nondimensional  Temperature  Excess.  1.2  Watts. 
Various  Spacings.  Top  Faces 
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Figure  42.    Nondimensional  Temperature  Excess.  1.2  Watts. 
Various  Spacings.  Bottom  Faces 
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Figure  43.    Nondimensional  Temperature  Excess.  1.2  Watts. 

Various  Spacings.  Heaters 
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Figure  44.    Nondimensional  Temperature  Excess.  0.2  Watts. 

Excl.  6.2  mm.  Front  Faces 
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Figure  45.    Nondimensional  Temperature  Excess.  0.2  Watts. 

Excl.  6.2  mm.  Right  Faces 


68 


CO 

□ 

6  X    i  -k3 

D 

□ 

O     |X<  + 

LEGEND 
=  0  MM 
=  14  MM 
=  10  MM 
=  23  MM 
=  NO  SHROUD 

CO 

□    O    <    +    X 

"  a 

q       ox  +   < 

-J 

> 

q         x  93- 

■* 

o 

!              n       ) 

CGW 

1                      o 

0   X-KJ 

G 

35 

_  d 

g; 

39 

92.* 

t                        i 

92*                        5iE 

92E                     92.2 

s-3  ♦i/na 

Figure  46.   Nondimensional  Temperature  Excess.  0.2  Watts. 
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Figure  47.    Nondimensional  Temperature  Excess.  0.2  Watts. 
Excl.  6.2  mm.  Bottom  Faces 
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Figure  48.    Nondimensional  Temperature  Excess.  0.2  Watts. 

Excl.  6.2  mm.  Heaters 
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Figure  49.    Nondimensional  Temperature  Excess.  1.2  Watts. 

Excl.  6.2  mm.  Front  Faces 
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Figure  50.    Nondimensional  Temperature  Excess.  1.2  Watts. 

Excl.  6.2  mm.  Right  Faces 


73 


Q2E 


tO< 


x      o    n  '<   + 


x     o       <   + 


n 

Q 

~  cc  3  2;  e  j  2: 
ii  ii  ii  ii  ii 

n  o  <  +  x 


xd     <  +'         a 


3  +0 


X     □     •« 


a      x  -ie! 


b     <ix 


00E 


91Z 

s-a  »i/na 


052 


-J 


522 


Figure  51.    Nondlmensional  Temperature  Excess.  1.2  Watts. 

Excl.  6.2  mm.  Top  Faces 
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Figure  52.    Nondimensional  Temperature  Excess.  1.2  Watts. 
Excl.  6.2  mm.  Bottom  Faces 
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Figure  53.    Nondimensional  Temperature  Excess.  1.2  Watts. 

Excl.  6.2  mm.  Heaters 
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yield  higher  temperatures  over  the  range  of  element  positions.  For 
higher  spacings,  there  was  no  consistency  in  the  relative  position  of 
the  curves  from  one  element  position  to  the  next.  This  indicates 
relative  independence  of  transport  on  channel  width  for  spacings 
greater  than  8  mm.  However,  data  for  1.2  watts  includes  even  the 
8  mm  spacing  within  the  mix  of  curves  (its  curve  does  not  consistently 
yield  higher  temperatures).  Therefore,  for  1.2  watts  the  indepen- 
dence of  transport  extends  below  the  8  mm  spacing.  General  patterns 
observed  for  both  powers  are: 

•  a  steady  rise  in  front  and  side  face  temperatures  with  element 
position  as  expected  due  to  exposure  to  increasing  convected 
energy  at  higher  positions. 

•  a  relative  levelling  of  top  and  bottom  temperatures  at  higher  posi- 
tions. This  is  perhaps  due  to  increasing  fluid  velocities  along  the 
channel. 

•  the  drop  in  heater  temperatures  at  block  7,  which  has  four  per- 
cent less  power  on  average. 

The  differences  between  the  two  powers  result  from  the 
stronger  driving  force  at  1.2  watts.  Increased  fluid  velocities  result  in 
improved  heat  transfer  characteristics,  making  8  mm  an  acceptable 
spacing.  This  point  is  emphasized  by  referring  to  the  graphs  of  com- 
ponent face  temperature  excess  vs.  shroud  spacing  for  0.2  watts 
(Figures  54-58)  and  1.2  watts  (Figures  59-63).  Significantly,  for  the 
1.2  watts  there  is  no  further  temperature  reduction  for  a  spacing 
greater  than  8  mm.  Indeed,  there  is  a  slight  rise  in  temperature  of 
the  latter  blocks  (7  and  8)  when  the  shroud  is  removed.    These  higher 
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Figure  54.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  0.2  Watts.  Front  Faces 
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Figure  55.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  0.2  Watts.  Right  Faces 
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Figure  56.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  0.2  Watts.  Top  Faces 
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Figure  57.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  0.2  Watts.  Bottom  Faces 
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Figure  58.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  0.2  Watts.  Heaters 
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Figure  59.   Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  1.2  Watts.  Front  Faces 
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Figure  60.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  1.2  Watts.  Right  Faces 
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Figure  61.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  1.2  Watts.  Top  Faces 
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Figure  62.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  1.2  Watts.  Bottom  Faces 
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Figure  63.    Nondimensional  Temperature  Excess  vs.  Shroud 
Position.  1.2  Watts.  Heaters 
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temperatures  (at  the  last  blocks)  are  not  observed  when  a  shroud  is  in 
position  because  the  global  forced  convection  within  the  channel  is 
driving  the  fluid  more  strongly  past  the  last  block.  As  seen  from  Fig- 
ures 54  through  58,  face  temperatures  at  a  power  level  of  0.2  watts  for 
8  mm  are  somewhat  higher  than  for  greater  spacings.  By  14  mm,  the 
difference  with  larger  spacings  begins  to  become  indistinguishable. 
Clearly  for  this  geometry  and  for  the  power  settings  examined,  a 
channel  width  of  2.5  times  component  depth  is  most  satisfactory 
while  a  width  of  1.3  times  element  depth  may  be  acceptable.  Further 
confirmation  of  the  observation  that  spacings  beyond  8  mm  do  not  sig- 
nificantly affect  heat  transfer  is  seen  in  Figures  64  through  68.  These 
show  the  Nu*  vs.  Gr*  data  along  with  the  correlation,  Nu*  =  1.88 
(Gr*)-15,  obtained  in  Chapter  III.  The  correlation  for  an  unshrouded 
surface  holds  very  well  through  14  mm.  At  8  mm  there  is  a  slightly 
greater  deviation.  As  anticipated,  it  is  inapplicable  for  the  smallest 
spacing  of  6.2  mm. 

B.     INDIVIDUAL  ELEMENTS  POWERED 

As  in  the  absence  of  a  shroud,  individual  blocks  were  powered  to 
1.0  watt,  one  at  a  time,  but  with  a  shroud  spaced  14  mm  from  the  sub- 
strate. The  resultant  data  is  displayed  in  Figures  69  through  73.  The 
plots  follow  similar  trends  as  for  the  no-shroud  condition  (Figures  21- 
25)  but  show  somewhat  smaller  temperatures.  Again,  the  protruding 
elements  have  minimal  effect  on  the  unheated  starting  length. 
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Figure  64.    Nu*  vs.  Gr*.  With  Correlation.  23  mm  Channel  Spacing 
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Figure  65.    Nu*  vs.  Gr*.  With  Correlation.  18  mm  Channel  Spacing 
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Figure  66.    Nu*  vs.  Gr*.  With  Correlation.  14  mm  Channel  Spacing 
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Figure  67.    Afo*  vs.  Gr*.  With  Correlation.  8  mm  Channel  Spacing 
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Figure  68.    Mi*  vs.  Gr*.  With  Correlation.  6.2  mm  Channel  Spacing 
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Figure  69.    Nondimenslonal  Temperature  Excess.  Front  Faces. 
Individually  Heated  Components.  14  mm  Channel  Spacing 
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Figure  70.   Nondimensional  Temperature  Excess.  Right  Faces. 
Individually  Heated  Components.  14  mm  Channel  Spacing 
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Figure  71.    Nondimensional  Temperature  Excess.  Tod  Faces. 
Individually  Heated  Components.  14  mm  Channel  Spacing 
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Figure  72.    Nondimensional  Temperature  Excess.  Bottom  Faces. 
Individually  Heated  Components.  14  mm  Channel  Spacing 
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Figure  73.    Nondimensional  Temperature  Excess.  Heaters. 
Individually  Heated  Components.  14  mm  Channel  Spacing 
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C.  COMBINATIONS  OF  ELEMENTS  HEATED 

In  a  manner  identical  to  that  without  a  shrouding  wall,  progres- 
sion of  elements  heated  was  made  up  the  array  while  all  upstream 
blocks  remained  powered.  Additionally,  block  8  was  powered 
throughout  and  a  shroud  14  mm  from  the  substrate  was  in  place.  The 
plots  (Figures  74-78)  reveal  that  for  this  channel  width,  as  with  no 
shroud,  the  powering  of  downstream  components  has  minimal  effect 
on  those  upstream  components  already  powered.  The  plots  further 
show  a  weak  effect  by  a  heated  element  (block  8)  on  those  compo- 
nents between  it  and  the  ones  powered  upstream,  which  are  heated 
only  by  the  convective  flow  generated  upstream  of  them. 

D.  ALL  ELEMENTS  POWERED,  TRANSIENT 
1.     Flow  Visualization 

Figure  79  shows  transient  flow  development  in  the  x-y  plane 
following  1.0  watt  power  step-up  with  a  19  mm  spaced  shroud  in 
place.  All  exposures  are  for  20  seconds.  The  photographs  from  left  to 
right  are  initiated  at  0,  20,  40,  100,  and  140  seconds  following  start  of 
heating.  In  the  first  picture,  parallel  flow  is  developing  over  blocks  3 
through  8  with  relatively  little  movement  at  blocks  1  and  2.  By  picture 
number  3,  the  flow  follows  more  closely  the  contours  of  the  protrud- 
ing blocks  and  originates  further  upstream.  By  photograph  5,  the  flow 
appears  nearly  fully  developed. 

The  power  step-down  from  0.2  watts  was  again  for  a  19  mm 
channel  width  and  20-second  exposures  (Figure  80).     The  pictures 
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Figure  74.    Nondimensional  Temperature  Excess.  Front  Faces. 
Multiple  Components  Heated.  14  mm  Channel  Spacing 
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Figure  75.    Nondimensional  Temperature  Excess.  Right  Faces. 
Multiple  Components  Heated.  14  mm  Channel  Spacing 
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Figure  76.    Nondimensional  Temperature  Excess.  Top  Faces. 
Multiple  Components  Heated.  14  mm  Channel  Spacing 
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Figure  77.    Nondimensional  Temperature  Excess.  Bottom  Faces, 
Multiple  Components  Heated.  14  mm  Channel  Spacing 
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Figure  78.    Nondimensional  Temperature  Excess.  Heaters.  Multiple 
Components  Heated.  14  mm  Channel  Spacing 
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Figure  79.    Flow  Response  During  1.0  Watt  Step-Up 
Transient.  19  mm  Channel 
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Figure  80.   Flow  Response  During  0.2  Watts  Step-Down 
Transient.  19  mm  Channel  Spacing 
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from  left  to  right  were  initiated  at  0,  60,  200,  320,  and  460  seconds, 
respectively,  following  cessation  of  power.  As  with  infinite  channel 
width,  the  flow  characteristic  persists  over  the  entire  period  observed. 
Only  in  the  last  picture  does  the  flow  begin  to  follow  less  closely  the 
contour  of  blocks  7  and  8  and  originate  farther  downstream. 
2.     Quantitative 

The  transient  temperature  response  curves  for  powers  of  1.0 
and  0.2  watts  and  a  shroud  spacing  of  19  mm  (Figures  81  and  82)  fol- 
low exactly  the  same  pattern  as  those  for  no  shroud  (Figures  12  and 
13).  For  the  spacing  examined,  the  shroud  does  not  affect  the  trans- 
port evolution. 
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Figure  8 1 .   Transient  Temperature  Response 
Following  Initiation  of  Heating.  1.0  Watt 
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Figure  82.   Transient  Temperature  Response 
Following  Cessation  of  Heating.  0.2  Watts 
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V.    CONCLUSIONS 

Flow  visualization  in  concert  with  quantitative  analysis  provided  a 
number  of  significant  determinations: 

•  The  convective  flow  was  three-dimensional. 

•  The  flow  was  laminar. 

•  The  thermal  layer  was  imbedded  within  the  momentum  exchange 
region. 

•  Heated    elements    have    minimal    effect    on    those    elements, 
upstream,  which  are  already  heated. 

•  There  is  no  gain  in  heat  transfer  for  channel  widths  greater  than 
two  times  component  depth. 

•  A  correlation  has  been  found  to  predict  heat  transfer  rates  for  the 
subject  geometry. 
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VI.    RECOMMENDATIONS 

It  is  suggested  that  the  following  areas  of  study  be  experimentally 
explored: 

•  Using  the  same  test  surface: 

-  explore  transient  response  for  other  block  faces 

-  with  the  traversing  thermocouple  probe,  determine  variation  of 
temperature  with  position  (in  the  thermal  layer)  for  all  blocks 

•  Construct  new  test  surface  and  rerun  the  same  battery  of  experi- 
ments to  determine  significant  changes  when: 

-  heated  elements  are  flush  mounted 

-  heated  elements  are  spaced  closer  in  the  x-direction. 
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APPENDIX  A 
SOFTWARE 

TEMPERATURE  ACQUISITION  PROGRAM 


10  M  II  I  I  I  I  I  I  I  I  I  I   I  M  II  I  I  II  M  I  |  |  I  I  I  |  I  |  I  I  |   I  |  |  I  I  |  M  I  I  M  I  I  I  I  M  I  I  |  I  I  I 

Z0  M  TEMPERATURE    ACQUISITION   PROGRAM  n 

30  i '  (STEADY    STATE)  ! ■ 

40  M  i  I  I  u  I  I  I  I  I  M  I  I  I  I  I  I  M  I  I  I  i  I  I  I  I  I  I  I  i  I  I  I  I  I  I  I  I  !  I  M  I  I  I  I  u  I  I  I  I  I  i  M  i 

50  REAL   Volts(60> 

G0  REAL    Te«p<59) 

70  CREATE    BOAT    "THES : .700,0, Z" ,48 

80  ASSIGN   6Path3    TO    "THES : ,700.0,2" 

90  PRINT    "  BLOCK    «8" 

100  PRINT    "    " 

110  OUTPUT   709i"CONFMEAS    DCV, 100, USE    0" 

1Z0  ENTER   709;Volts<60) 

130  OUTPUT   703i"CONFMEAS   DCV . 100- 105 .USE    0" 

140  FOR    1-0   TO    5 

150  ENTER   709jVolts< I  ) 

160  Te*p( I )-. 0006797+ (Z58Z 5. 13Z8»VoltsU  > >-< 607789 . Z467»< Vol ts< I >»Volts< I ) ) )-( 

Z 195Z034. 3364* < Volt s<  I  >"3>  >  +  (  8370810996.  1  874* (  Volt s<  I  )*4)  ) 

170  PRINT    "T.    C.       r'sl  +  l,"       Volts    D.C".       "  !  Volts<  I  )  .  "Tenp.    DEG.    C    "  ;  Tenp<  I  ) 

180  NEXT    I 

190  PRINT    "    " 

Z00  PRINT    "  BLOCK    tt7" 

Z10  PRINT    "    " 

ZZ0  OUTPUT    709;"CONFMEAS    DCV . 106- 1 1 1 . USE    0" 

Z30  FOR    I-G    TO    11 

Z40  ENTER   709;Volts< I ) 

Z50  Te«p< I )-. 0006797  +  < Z58Z 5. 13Z8«Volts(  I ) >-< 607789 . Z467»< Vol ts< I ) •Volts( I ) ) )-( 

Z 195Z034. 3364* < Volt s< I)*3) >  +  ( 83708 1099G. 1 874* < Volt s< I  )"4)  ) 

ZG0  PRINT    "T.    C.       »";I+1,"      Volts    D.C.       " ? Vol ts< I  )  , "Tewp.    DEG.    C    " ; Tenp( I ) 

Z70  NEXT    I 

Z80  PRINT    "    " 

Z90  PRINT    "  BLOCK    »6" 

300  PRINT    "    " 

310  OUTPUT  709i"CONFMEAS  DCV , 1 1 Z- 1 17 .USE  0" 

3Z0  FOR  I-1Z  TO  17 

330  ENTER   7095Volts< I) 

340  Tenp< I )-.0006797  +  <  Z58Z5.  I 3Z8»Vol ta(  I ) >-< 607789 . Z467*< Vol ts< I )• Vol ts( I ) ) )-( 

Z19SZ034.3364»<Volts< I  )*3) )+( 8370810996. 1 874* ( Volt s( I )*4) > 

350  PRINT    "T.    C.       rtl  +  l,"      Volts    D.C.       "  5  Volts(  I  )  .  "Temp.    DEG.    C    "  i  Temp<  I  ) 

360  NEXT    I 

370  PRINT    "    " 

380  PRINT   "  BLOCK    »5" 

390  PRINT    "    " 

400  OUTPUT  709;"CONFMEAS  DCV , 1 1 8- 1 1 9 , USE  0" 

410  FOR  1=18  TO  19 

4Z0  ENTER   709;Volts( I ) 
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4  30  Temp(  I  ) -  .  0006797+ (  Z58ZS .  13Z8»Volts<  I  )  )-  <  607789  .  2467* (  Vol  ts<  I  )  •Vol  ts<  I  )  )  )-( 

Z19SZ034. 3364«< Volts( I >"3)  >  +  < 83708 10996. 1 874»( Volts( I  )"4)  ) 

440  PRINT    "T.    C.       r';I+1,"      Volts   D.C.       " ; Volis( I ) ,"Tenp.    DEG.    C    " ; Tenp< I ) 

450  NEXT    I 

4G0.  OUTPUT  709i"CONFMEAS  DCV , Z00-Z03 , USE  0" 

470  FOR  I-Z0  TO  Z3 

480  ENTER  709;Volts<  I  > 

490  Temp< I >■ . 0006797+< Z58Z5. !3Z8*Volts(  I ) )-( G07789 . Z467M Vol ts< I )»Volts( !)))-< 

Z195Z034. 3364* < Vol ts( I r 3 >>+< 8370810996. 1874MVolts< I)*4)) 

500  PRINT  "T.  C.   »"?I+1,"   Volts  D.C.   "  s  Vol  ts(  I  )  ,  "Tetnp.  DEG.  C  "  i  Tenp(  I  > 

510  NEXT  I 

5Z0  PRINT  "  " 

530  PRINT  "  BLOCK  *4" 

540  PRINT  "  " 

550  OUTPUT  709?"CONFMEAS  DCV .Z04-Z09 ,USE  0" 

560  FOR  I-Z4  TO  29 

570  ENTER  709;Volts(I) 

580  Tenp( I >- .0006797+< Z5825. 13Z8«Volts(I ) >-< 607789 . Z467»< Volt s< I ) •Volts* I ) )  )-( 

Z195Z034.3364«< Volts( I )* 3) )  +  (  8370810996. 1874* < Volts( I )A4) ) 

590  PRINT  "T.  C.   *"jl+1,"   Volts  D.C.   " ? Vol ts(  I ) , "Temp.  DEG.  C  ";Tomp(I) 

600  NEXT  I 

B10  PRINT  "  " 

6Z0  PRINT  "  BLOCK  *3" 

630  PRINT  "  " 

640  OUTPUT  709;"CONFMEAS  DCV , Z 10-Z 1 5 .USE  0" 

650  FOR  1-30  TO  35 

660  ENTER  709»Volts(I) 

670  Tenp< I )= .0006797+< Z58Z5. 13Z8«Volts<  I ) )-< 607789 . Z467»< Vol ts< I )»Vol ts< I ) ))-( 

Z195Z034. 3364* (  Vol ts(  I )A3))  +  (  837081 0996. 1874»<Volts( I >A4>> 

B80  PRINT  "T.  C.   K"»I+1,"   Volts  D.C.   " ; Volts( I ) , "Temp.  DEG.  C  "iTenp(I) 

690  NEXT  I 

700  FOR  J- J  TO  14  - 

710  PRINT  "  " 

7Z0  NEXT  J 

730  PRINT  "  BLOCK  ttZ" 

740  PRINT  "  " 

750  OUTPUT   709j"CONFMEAS   DCV , Z 1 6-Z 1 9 ,USE    0" 

760  FOR   1-36   TO    39 

770  ENTER  709;Volts<I) 

780  Tenp( I )- .0006797+ (  Z58Z 5. 13Z8»Volts(  I ) )-< 607789 . Z467M Vol ta< I MVolts< I ) )  )-< 

Z195Z034. 3364* < Vol ts(  I >" 3) )+< 8370810996. 1874»(Volts( I )*4)> 

790  PRINT  "T.  C.   r'|I  +  1,"   Volts  D.C.   " ; Vol ts< I ) . "Tenp.  DEG.  C  " ; Tenp< I  ) 

800  NEXT  I 

810  OUTPUT  709i"CONFMEflS  DCV , 300- 301 , USE  0" 

8Z0  FOR  1-40  TO  41 

830  ENTER  709;Volts(I) 

840  Te«p< I >  =  .0006797+(Z58Z5. 13Z8»Volts(  I ))-( 607789 . Z467»(  Volt s< I )»Vol ts( I ) ) )-( 

Z195Z034.3364»(Volts( I >" 3 )>  +  (  8370810996. 1874»<Volts( I )*4)> 

850  PRINT  "T.  C.   «"il+1."   Volts  D.C.   " ; Volts< I ) . "Tenp.  DEG.  C  "jTemp<I) 

860  NEXT  I 

870  PRINT  "  " 

880  PRINT  "  BLOCK  *\" 

890  PRINT  "  " 

900  OUTPUT  709i"CONFMEAS  DCV . 30Z-307.USE  0" 

910  FOR  I-4Z  TO  47 

9Z0  ENTER  709;Volts<I) 
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9  30        Tenp( I )=.0006797+<  25825. 1 32  8»Vol ts< I ) )-( 607789 . 2 467»( Vol ts(  I ) •Vol ts( I ) ) >-( 
Zl 952034. 3364»< Vol ts< I ) A 3 ) ) +< 8370810996.  1874M Volts< I  )*4)  ) 

T.    C.       r'iI+1,"      Volts    D.C.       "iVolts< I > ."Temp.    DEG.    C    " ; Temp( I > 


ePath3;Tenp( • ) 


BATH   TEMPERATURES    (TOP    TO   BOTTOM)" 

709;"CONFMEAS   DCV , 317-31 9 .USE   0" 
57   TO   59 
>9iVolts<  I  ) 

.0006797+< 25825. 1328»Volts(  I ) )-< 607789 . 2467* ( Vol ts< I >*Volts< I > > )-< 
2 1952034. 3364* (  Volt s< I )A 3) )+< 8370810996. 1874»(  Volts<  I >*4>) 
1090      PRINT    "VOLTS    D.C.       " ; Vol ts( I > , "TEMP .    DEG.    C    " ; Temp( I > 
1 100      NEXT    I 
1110      ENO 


940 

PRINT 

950 

NEXT  I 

9  60 

! 

970 

I 

9  80 

OUTPUT 

990 

I 

1000 

! 

1010 

PRINT 

1020 

PRINT 

1030 

PRINT 

1040 

PRINT 

1050 

OUTPUT 

1060 

FOR  I-1 

1070 

ENTER  ' 

1080 

Te«p< I 
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DATA  REDUCTION  PROGRAM 


10  M  I  I  II  I   I  I  I  I   I  I   I  M   II   I   I   I  M  I   M  I  I  I  I  I  I  !  I  I  I  I   I   I  |  |  I  I  I  I  I  I  I  I  I  I  I  I  M  |  I  M  I  I   I   I   I  I  I 

Z0  !  I  •  ! 

30  i  i  DATA   REDUCING   PROGRAM    (BASIC)  !! 

40  i i  !  i 

50  I  I  I  I  I   I  I  I  I  I  I  I  M   I  I  II  I   I  I  I  II  I  I  I  I  I  I  I  t  I  I  I  I  I  I  I  J  U  I  I   I   I  I  I  I   I  I  I  I  M  I  I  I  I  I  I   I   I  I   I   I 

G0  REAL    Tewp(48> ,  Tt<48) ,Qcond(8) ,Qconv(  8) ,Qf lux< 8 ) , H( 8 ) .Nu(8) ,Nustar<8) 

70  REAL   Gr(8)  ,Ttt(8)  ,  At<  8)  .Dtf  (  8)  .Dtr<  8  )  ,Dil(  8  >  ,Dit(8)  ,Dtb<8)  ,Dth<8)  ,B1(8) 

80  REAL    Jstar(8),Jp(8)  ,Gratar(8) 

90  REAL   Pouer .Space, Dxpg.Dxr.Kpg.Kr .Block 

100  REAL    A.G.Dis.Try , KhZo .Beta.Kinvisc ,Pr 

I  10  INTEGER    I.J.K.L.M.N.O 

1Z0  ! 

130  Dxpg«. 006731 

140  Dxr- .003175 

150  Kpg-.14Z1 

160  Kr- .0389 

170  Chpht-. 00787 

180  Chpdpth-.00610 

190  Chpudth-.0Z388 

Z00  PriR-(Z»(Chpht+Chpudth))+(Z»(Z»(Chpht+Chpdpth)))+(Z»(Z»(Chpudth+Chpdpth))) 

Z10  A1«Chpudth»Chpht 

ZZ0  AZ=Chpudth»Chpdp-th 

Z30  A3=»Chpht»Chpdpth 

Z40  Atotal»Al+<Z.»AZ)+(Z.»A3> 

Z50  Atotala-(Z.*AZ)+(Z.»A3) 

Z60  G=9.81 

Z70  M   II  M  I   I  I  I   |  I  I   I   I  I  I   I   I   I   I  I  I  I  I  I  I   I  I  M  I   I  M  I   I   I   I   I   I  I  !   I   I   I   I  I  |  I   I   I  I  I  I   I   I   I  I   I   I   M   I   I   I   I  I   I   I   I   | 

Z80  '!  CALCULATE    CHARACTERISTIC    LENGTH  !! 

290  I  I  I  I i  I  I  !  I |  i  i  l  I  l  !  I  t  I  i  i  I  |  I M  i  i I  I  I l  i  mi  i  |  I  |  i  I  M  I  Mi  I  I  l  l  I  I  I  |  I  i  I  i  M  I  I  I  I  !  I  |  !  M 

300  Lbar=A-to-tal/Prim 

310  I 

3Z0  R1=Dxpg/(Kpg»Al  ) 

330  RZ»Dxr/(Kr»A1 ) 

340  I 

350  PRINTER  IS  1 

360  PRINT  "INPUT  RUN  8,  INTEGER  ONLY" 

370  INPUT  Try 

380  PRINT  "  " 

390  PRINT  "  " 

400  PRINT  "INPUT  POUER  VALUE  IN  UATTS" 

410  INPUT  Power 

4Z0  PRINT  "  " 

430  PRINT  "  " 

440  PRINT  "INPUT  BLOCK  8'S  POUERED" 

450  INPUT  Blnp$ 

460  PRINT  "  " 

470  PRINT  "  " 

480  PRINT  "INPUT  SPACING  IN  MM" 

490  INPUT  Spaces 

500  PRINT  "  " 
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510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
7Z0 
730 
740 
750 
760 
770 
780 
790 
800 
810 
8Z0 


PRINT  "  " 

PRINT  "INPUT  AMBIENT  TEMPERATURE" 

INPUT  Tinf 

PRINT  "  " 

PRINT  "  " 

i   PRINT  "INPUT  NEU  FILENAMES" 

l   INPUT  NeufileS 

!   PRINT  "  " 

!   PRINT  "  " 

PRINT  "INPUT  OLD  FILENAMES" 

INPUT  OldfileS 


!  CREATE  BDAT  NeufileS, 128 
i  ASSI6N  SRoad  TO  NeufileS 
I 

ASSIGN  flStreet  TO  OldfileS 

ENTER  eS-treet;Tenp(0 

FOR  1-0  TO  47 

J-J+1 

! PRINT    J,Tenp< I  ) 

NEXT  I 


M  I  I  M  M  M  I  I  I  I  I  I  I  M  M  I  I  I  I  I  |  I  I  I  I  II  |  |  |  I  I  I  I  I  I  I  I  II  I  I  |  |  I  I  M  I  I  I  I  It  I  I  I  I  I 

! i  CALCULATE  AVERAGE  SURFACE  TEMPERATURES  ! ! 

I  |  |  |  |  I  |  |  i  I  I  i  |  |  |  |  !  I  I  I  I  I  i  I  i  |  I  I  i  |  i  |  I  I  i  |  I  I  I  i  M  I  I  i  I  j  j  i  M  i  i  ii  I  ii  |  m  I  I  |  i 

•    PRINT  "AVG.  SURFACE  TEMPERATURE" 
J-0 

FOR  1-1  TO  8 
IF  1-6  THEN 

Ai( I )  =  ( ( Temp< J+1  )»A3>+< Tenp<J+2 )»A3)  +  <  Temp< J+3)»A2 )+<Tenp<  J+4)»A2) >/ 
Atotala 
830  ELSE 

840  At( I )-< ( Temp< J)»A1  >  +  ( Tenp<  J+l >»A3>+< Tenp< J+2 )»A3)+(Temp<  J+3)«A2)+( Te 

mp<  J+4)»A2>  >/A+.o+.al 


850 
860 
870 
880 
890 
900 
910 
920 
9  30 
940 
950 
960 
970 
980 
990 
1000 


END  IF 
i   PRINT  A-t(  I  > 
J-J+6 
NEXT  I 
I 
i  i  i  i  i  it  |  i  i  ii  i  i  I  i  i  i  I  i  i  i  i  ii  I  i  I  i  i  i  I  i  i  i  I  i  i  i  i  |  i  i  i  i  i  i  i  i  ii  i  i  i  i  i  i  i  i  I  i  i  i  i  i 

i  i  CALCULATE  TFILM  ■  ! 

i  i  ii  i  i  i  i  ii  i  i  i  i  I  i  i i i i  i  i ii  i  i  i  i  |  i  i  i  I  i  |  I  I  i  I  |  I i  i  i  i  i  i  i  i  i  i  i  i  i  i  ii  i  |  i  I  i  i  it 

J-1 

Tbar-< At( J)+At( J+1  ) +At (  J  +  2  ) +At (  J+3)+At(  J+4  )+At<  J  +  5  )+A+.(  J  +  6)+At  (  J  +  7)  )/8.0 

Tf ilnc=< ( Tbar+Tinf )/2.0) 

Tf iln-Tf ilnc+273. 15 

PRINT    "TFILM   DEGREES   C" ,Tf ilnc 

PRINT   "TFILM". Tfiln 
l  i  |  I  I  I  I  I  I  M  i  I  I  I  I  I  I  I  I  I  I  I  i  i  I  I  I  I  I  I  I  I  |  I  |  I  (  I  i  11  I  I  I  j  i  I  I  I  i  I  I  I  I  I  I  I  I  I  I  I  I  I 

!  CALCULATE    UATER  PROPERTIES  !  ! 
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1010  i  i  i  i  i  i  i  i  i  i  i  I  I  !  I  I  M  I  M  M  !  I  I  !  M  I  I  I  I  I  I  !  I  I  I  I  I  I  I  I  I  I  M  M  M  I  I  I  I  I  I  I  !  I  I  I  M 

10Z0  8eta-< -7.9448548E-8»<Tf iIpTZ) >+< 5. 7356479E-5«Tf i In)- .00978 10563 

103®  Spvol-(4.69699E-9»(Tf iln*Z>)+< -Z . S3745E-6«Tf l In) +  . 001341 903 

1040  Dynvisc-<  3. Z 34851 1E-7M Tf iln'Z ) )  +  ( - .000Z 1 474487* Tf iln)*  .0361 G679Z 

1050  K invisc-Spvol *Dynvisc 

1060  KhZo-<  1  .4l8181818E-3»Tf iln)+. 186G 

1070  Pr-<4.65706E-3»(Tf ilm*2>>+< -Z . 9ZZ094»Tf iln >+463. 331 9 

1080  M  I  I  I  I  i  M I  i  I  I  i  I  I  i  I  I  I  M  i  i  I  I  i  I  I  I  |  I  M  I  I  I  M  i  M  I  I  M  I ! i  |  I  I  I  I  M  I  i  !  I  |  I  I  I  I 

1090  !!  CALCULATE  TEMP-TINF  M 

1 100  I i  i I  I  I  I  I  M  I  I  I  !  I  I  I  i  M  I  i  I  I  ! I  I  i I  I  I  I  !  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  i  I  I  I  |  I  I  |  I  I  I  I  I  I  I  I  I 

1  1  10  !    PRINT  "TEMP-TINF" 

1  1Z0  FOR  1=0  TO  47 

11 30     IF  1-5  THEN 

1140       Tt(  D-888.888 

1150     ELSE 

1160        IF  1-30  THEN 

1170  Tt(  D-888.888 

1180        ELSE 

1190       Tt( I)-Tenp< I)-Tinf 

1Z00     END  IF 

1Z10     END  IF 

1ZZ0  i   PRINT  Tt< I ) 

1Z30  NEXT  I 

1 Z  40  i  i  i  I  I  I  i  I  I  I  I  I  I  |  I  I  I  I  |  I  I  I  I  M  I  i  i  I  !  I  M  l  M  M  I  I  !  |  I  I  I  |  I  I  I |  I  I  I  | i  I  I  I  M  M  |  I  I  I 

1Z50  !!     CALCULATE  CONDUCTION  LOSSES  THROUGH  THE  TEST  SURFACE       !! 

1  Z  60  I  i  i  i  I  i  I i  I  I  I  I  I  i|  I  I  I  I  I  M  I  I  I  I  I  I  I  I  I  I  I  i  M  M  I  i  I  I  I  M  I  I  I  I  I  M  i | I  I  I  I  |  I  I  I  I  i  i  I 

1Z70  !    PRINT  "QCOND" 

1Z80  L-5 

1Z90  FOR  1=1  TO  8 

1300     IF  1=1  THEN 

1310         Qcond( I )-.015»Pouer  *- 

13Z0     ELSE 

1 330        Qcond( I  )-Tt( L  )•(  1 .0/< Rl +RZ ) ) 

1340     END  IF 

1350  i   PRINT  Ocond( I) 

1360  L=L+6 

1370  NEXT  I 

1  380  I  I i i  I  I  I  i  I  I  I  I  !  i  i  I  I  I  i  I  i  I  i  I  I  i  u  u  I  I  i  u  i  I  u  I  i  i  i  t  I  i  I  I  I  i  I  i  I  i  i  i  I  I  u  I  i  i  I  i  I  i 

1390  it  CALCULATE  CONVECTED  HEAT  FLUX  !! 

1 400  i  I  I  I  I  I  I  I  |  I  I  I  I  I  I  |  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  M  I  I  M  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  i  I  I  M 

1410  i    PRINT  "QCONV 

14Z0  FOR  1-1  TO  8 

1430  Qconv< I )=Pouer-Qcond< I  ) 

1440  Qf lux(  I )-Qconv< I >/A1 

1450  NEXT  I 

1 4G0  i  I  I  I  i  ii  i  I  i  |  i  I  i  I  i  i  I  I  i  i  i  i  i  I  i  i  i  I  i  i  i  I  I  I  i  i  i  I  I  |  I  i  |  I  i  i  i  I  I  i  I  I  |  i  I  I  i  i  i  i  i  I  i  I  i  i 

1470  M  CALCULATE  TEMPERATURE  SCALING  FACTOR  II 

1 480  I  i ! I  i  i  i  I  I  i  !  i  I  l  I  I  I  i  i  i  i  i  I  I  I  I  i I  I  l  I  I  I |  I  I  I  I  i  I  I |  I  I  I i I  I  i  i  I  i I  I  M  i  I  i  i  I  I  I  i  I  i  I 

1490  I   PRINT  "J*" 

1500  FOR  1-1  TO  8 
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1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1630 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
19  30 
1940 
1950 
1960 
1970 
1980 
1990 
2000 


TstaH I )-Qf lux( I )»Lbar/Kh2o 
!    PRINT  TSTAR( I) 
NEXT  I 

I  I  I  M  II  II  II  I  I  I  I  I  I  I  II  I  I  I  I  I  I  II  I  !  I  I  I  I  I  I  M  II  I  II  I  I  II  I  I  I  I  M  I  I  I  II  M  I  I  M  I  I 

!  !  CALCULATE    MODIFIED  GRASHOF   NUMBER  ! 

I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  |  I  I  I  I  I  I  I  I  I  I  I  I  II  I  I  I  I  I  I 

FOR  1-1  TO  8 

Srstar< I )-G«Beta*Qf lux( I )•< Lbar*4)/(  Kh2o»< K invisc*2 > ) 
NEXT  I 
I  I  I  I  II  I  I  II  I  I  I  I  I  I  I  I  I  (  I  II  I  |  I  I  I  I  I  I  I  I  I  II  I  I  I  I  |  I  M  I  (  I  II  M  I  |  I  I  I  I  I  I  I  I  I  I  I  I  |  I 

i  !  CALCULATE  HEAT  TRANSFER  COEFFICIENT  ! ! 

I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  i  I  I  I  I  I  I  I  I  I  I  i  i  I  I  I  I  I  I  I  I  I  I  I  !  I  I  I  I  I  I  I  i  I  i  I  I  I  I  I  j  I  I 

FOR  1-1  TO  8 

H(I >-Qflux( I )/(At( I )-Tinf ) 
NEXT  I 

I  I  I  I  I  I  I  I  I  I  I  I  |  I  I  I  !  M  I  I  I  I  !  I  I  !  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  ;  I 

!  CALCULATE   NUSSELT   NUMBER  !  ! 

M  l  I  M  I  I  |  |  |  I  M  i  |  I  I  I  I  I  j  i  I  I  I  i  I  i  i  I  i  i  i  i  i  I  I  i  i  j  I  I  I  |  l  !  I  M  I  !  I  I  I  M  I  I  I  i  I  I  i  |  |  I 

FOR    1-1    TO    8 

Nu( I )-H(I )»Lbar/Kh2o 
NEXT  I 
M  i  i  i  i  i i i |  i  | i  i i i  i  i  i  i  |  i  i  i  |  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i i i i  i  i  i  I  |  i  i  i  i  i  i  i  i  i  i  i  i  i  i  |  ii  i  i  i 

M  CALCULATE  NUSSELT  SCALING  FACTOR  II 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  |  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  (  i  |  i  i  i  i  |  |  i  i  i  i  i  i  i  i  i  i  i  | 

FOR  I-l  TO  8 

Js-9-I 

Nus-tar(  I  )-Nu(  I  )»Js*(  1/6) 
NEXT  I 

I  |  |  |  I  II  |  |  I  |  |  I  I  i  j  |  I  I  I  I  I  |  |  I  I  I  I  i  I  I  I  I  I  I  M  I  i  I  i  |  I  I  |  l  |  I  ||  |  !  |  |  i  I  i  i  i  I  I  M  j 

!  GENERATE  OUTPUT  DATA  FILES 

i  i  i  l  i  |  i  i  |  i  ii  I  |  M  |  I  |  I  I  I  i  i  i  i  i  I  i  i  i  i i i  i  |  |  i|  |  (  I  I  i|  I  i  i  |  i  i  i  |  i  i  i i  i  i  ||  i  i  i 

PRINTER    IS    701 
I 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  "  " 

PRINT  "  " 

PRINT  "  " 

PRINT  "BLOCK  »" , "QCOND" , "QCONV" , "GR  »" ,"T»" ,"H" ,"NU  »"."NU»" 

FOR  1=1  TO  8 


'"RUN    NUMBER"" , 15X. DO";  Try 
'"POUER   IN   UATTS"" , 1 1 X , Z. D" iPouer 
'"SPACING    IN   MM"" , 1ZX, 10A"; SpaceS 
'"BLOCK    tt's   POUERED"",8X. 10A";Blnp$ 
"AMBIENT   TEMP.    DEG.    C"" ,5X ,DD. DO" i Tinf 
'"FILM    TEMPERATURE "", 9 X.DDO. DO" ;Tf i In 
'"FILM    TEMPERATURE,    C"" ,6X ,DD.  DD" ; Tf l lnc 
'"THERMAL    COND.    OF   UATER" " , 3X , Z . DOD" i KhZo 
'"EXPANSION   COEFFICIENT^"". 2X,3DE"»Beta 
'"KINEMATIC    VISCOSITY"" ,6X , 3DE" ; K invisc 
'"PRANDTL    NUMBER"" , 1 lX,DD.DD";Pr 
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Z010 

2020 

DE" tBl 

2030 

2040 

2050 

2060 

2070 

2060 

2090 

2100 

21  10 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

Dif(I) 

2250 

22G0 

2270 

2280 

2290 

2  300 

2310 

2320 

2  330 

2  340 

2350 

2  3G0 

2  370 

2  380 

2  390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 


Block-9-I 

PRINT  USING  " 3X , D, 3X , S3DE , 3X , SZ . DOO , 3X . S3DE . 2X , S3DE , 3X , S3DE , 2X , S3DE , 2X . S3 
ock ,0cond(  I ) ,0conv(  I  ) ,6rstar<  I  ) ,Tstar<  I ) ,H< I ) ,Nu< I ) ,Nustar(  I ) 
NEXT  I 
PRINT  "  " 
PRINT  "  " 
PRINT  "  " 
PRINT  "  " 

I   I  M   I  I  I  !  I  M   I   I  I  I  I  M  I  I  I  I  II   |  I  I  M  II   I  M  I   I   M   I  |  II  II  I  I  I  II  I  I  I  I  I  I  I  II  II   I   I  II  I  I  I  I  I 

I!  CALCULATE    TEMPERATURE    DIFFERENCES  !! 

I  M  I  I  I  I  i  I  M  '  •  I  I  I  i  I  I  I  !  I  I  I  I  I  !  i  !  M  I  I  I  I  !  I  M  I  I  I  I  I  I  I  M  |  !  I  !  I  I  |  I  I  |  i  i  |  i  i  i  |  |  i  i  i 

K-0 

FOR   N-1    TO   8 

Dif(N)-Tt(K) 

Dtr(N>-Tt(K+1 ) 

Dtl(N)-Tt(K+Z) 

Dtt(N)»Tt(K+3) 

Dtb(N)-Tt(K+4) 

Dth(N)-Tt(K+5) 

K-K+6 
NEXT   N 

PRINT    "BLOCK    «" ,"OTF" , "DTR" , "DTL" , "DTT" , "DTB" , "DTH" 
FOR   1-1    TO   8 
Block-9-I 

PRI NT   USI NG   " 3X , D , 3X . S3DE . 3X , S3DE , 3X , S3DE . 3X , S3DE . 2 X , S3DE , 2X , S3DE" ; Bl ock , 
,Dtr(I>,Dtl< I).Ott(I) ,Dtb< I),Dth(  I ) 
NEXT    I 
l  u  u  i  I  I  I  M  i  I  I  I  I  I  i  I  I  I  l  I  l  I  I  H  !  i  u  u  I  I  i  I  i  I  I  I  M  i  I  I  I  i  I  I  I  i  I  u  I  u  I  i  I  I  I  I  I  I 

!!  CALCULATE   NON   DIMENSIONAL    TEMPERATURE    DIFFERENCES  ! 

I  I  j  j  M  I  i  I  i  |  I  I  !  !  |  I  (  I  |  I  t  I  j  I  !  j  I  I  I  |  M  I  (  i)  i  I  I  I  I  I  |  I  I  |  i  i  M  i  I  I  I  i  |  i  I  i  i  I  I  I  I  i  i 

FOR    1=1    TO    8 

i  IF    1-1    THEN 

I  Dthn( I )-88888 

i  ELSE 

i         IF    1-6    THEN 

!  Dtfn(  I )-88888 

Cttfn<  I)=D-tf<  I  )/Tstar<  I 
Dtrn(  I )-Dtr< I)/Tstar( I 
Dtln<  I>*Dtl<I)/Tstar<I 
D-ttn(  I>-Dt"t<  I)/Ts-tar(  I 
Dtbn<  I  )-D-tb<  I  )/Ts-tar(  I 
Dthn(  I  )=Dth(  I  )/Tatar(  I 
IF    1-1    THEN 
D-thn(  I  )-88888 
ELSE 

IF    I-G   THEN 
Dtfn< I >-88888 
END    IF 
END   IF 
NEXT    I 
PRINT   "    " 
PRINT    "    " 
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2510  PRINT    "X/L"  ."DTFN" ."DTRN" ,"OTLN" ."OTTN" ,"DTBN" ,"OTHN" 

2520  Xs-.795 

2530  6ap-.0ZS4/.3048 

2540  FOR   1-1    TO   8 

2550  Xl-Xs-<(I-1 )»6ap) 

2560  PRINT   USIN6    " . DOO, 3X ,S3D€ , 3X , S3DE . 3X . S3DE . 3X , S3DE ,2X , S3DE , 2X ,S3DE" ; XI ,Dtf 

n(  I  )  .Dtrn(  I  >  ,D-tln(  I  >  .Dttn<  I  )  ,0"tbn(  I  )  ,Dthn<  I  ) 

2570  NEXT    I 

2580  PRINT    "    " 

2590  PRINT    "    " 

2600  PRINT    "DTR   -    RI6HT    T.    C.    -    TINF.    etc." 

2610  PRINT   "T«    IS    TEMPERATURE    SCALING   FACTOR" 

26Z0  PRINT   "GR   »    IS   MODIFIED  GRASHOF    NUMBER" 

2630  PRINT    "H   IS   HEAT    TRANSFER   COEFFICIENT" 

2640  PRINT    "NU   «    IS    NUSSELT    NUMBER" 

2650  !  OUTPUT   BRoadiQcond< * ) ,Qconv<* ) ,6r< • ) ,Ttt< * > ,B1( * > ,Dtf ( • ) ,Dtr< • ) , 

Dtl<»).Dtt<0  .Dtb(*)  .Dth(  •) 

2660  i 

2670  I i   NOTE  THIS  PROGRAM  IS  USED  IN  CONJUNCTION  UITH  THE 

2680  !!   THE  DATA  AQUISITION  PROGRAM  AND  ACCOUNTS  FOR  THE 

2690  I!   INOPERATIVE  THERMOCOUPLES  ON  BLOCKS  3H8 

2700  i 

2710  END 
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THERMOCOUPLE  VARIATION  WITH  TIME  PROGRAM 


10  M  l|  II  |  I   I  I  I  I   I   I  I  I  I  II   I  I  I  |  I  I   I  |  |  I  I   II  II  I  I  M  II  I  M  I  I  I  I  |  I   M  I  I  I  I  II   I  II   I  I   J   I 

20  i I  TEMPERATURE    TRANSIENT   RESPONSE   PROGRAM  ! I 

30  i  i i  i i i  i  |  i  i  i  i  i  I  i  M  M  M  M  i  i  I  |  I  I  I  I  I  I  i i  !  M  I  I  I  I  I  I  I  I  I  I  I  I  i  i  i  i  i  i  i  i  m  i  i  i  I 

40  REAL    TZ(575),T3(575) , T4( 575 > , T5< 575 ) . T6< 575) ,T7< 57S > , T8< 575 > 

50  REAL   VZ(S75>  ,V3<575>  ,V4<575)  ,V5<575)  ,V6<575),V7<575)  ,V8(57S> 

60  REAL   G(575) 

70  N-575 

80  PRINTER  IS  701 

90  CREATE  BOAT  "TRANS : ,700,0, 1" ,502 

100  ASSIGN  ©Trail  1  TO  "TRANS :, 700. 0,1" 

110  II  I  I  II  I  I  M  I  I  J  I  I  I  I  I  I  I  I  I  |  I  I  I  I  II  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  II  I  I  I  I  I  I  I  I  I  I  I  M  I  I  I  I  I  I  I  II  I  I  I  I  I 

120  BEEP 

130  FOR  J-l  TO  5 

140       OUTPUT  709;"CONFMEAS  DCV USE  0" 

150        ENTER  709s VI (J) 

160      OUTPUT  709j"CONFMEAS  DCV USE  0" 

170        ENTER  709}V2(J) 

180      OUTPUT  709;"CONFMEAS  DCV USE  0" 

190        ENTER  709;V3(J) 

200  OUTPUT   709;"CONFMEAS   DCV USE   0" 

2  10  ENTER   709;V4(J) 

220  OUTPUT   709;"CONFMEAS   DCV USE   0" 

2  30  ENTER   7095V5(J) 

240  OUTPUT   709j"CONFMEAS   DCV USE   0" 

250  ENTER  709;V6(J) 

260  OUTPUT   709;"CONFMEAS   DCV ,USE   0" 

270  ENTER  709;V7<J)  *- 

280      OUTPUT  709;"CONFMEAS  DCV USE  0" 

290        ENTER  709)V8<J) 

300        NEXT  J 

310        BEEP 

320        FOR  J-6  TO  N 

330      OUTPUT  709;"CONFMEAS  DCV USE  0" 

340         ENTER  709; VI (J) 

350  OUTPUT   709?"CONFMEAS   DCV USE   0" 

360  ENTER   709;V2(J) 

370  OUTPUT   709;"CONFMEAS   DCV USE   0" 

380  ENTER   709;V3(J) 

390  OUTPUT   709;"CONFMEAS   DCV USE    0" 

400  ENTER   709-,V4(  J) 

410      OUTPUT  709j"CONFMEAS  DCV USE  0" 

420         ENTER  709;V5< J) 

430       OUTPUT  709;"CONFMEAS  DCV USE  0" 

440         ENTER  709;V6< J) 

450      OUTPUT  7095"CONFMEAS  DCV USE  0" 

460         ENTER  709;V7( J) 

470       OUTPUT  709;"CONFMEAS  DCV USE  0" 

480         ENTER  709;V8( J) 

490         IF  J-N  THEN 

500        BEEP  2450.97»Tone,1 .5 

510        END  IF 

520        NEXT  J 
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530  A=. 0006797 

540  B=Z58ZS.13Z8 

550  C--607789.Z467 

560  D--Z195Z034. 3364 

570  E=8370810996. 1874 

580  FOR  J-1  TO  N 

590  IF  J<5  THEN 

600  G(J)-0 

610  ELSE 

6Z0  G( J)-< J-5)»Tim 

630  I   Tin  IS  TIhE  BASED  ON  TIME  PROGRAM  RUNS  FROM  ZND  BEEP  TO  TONE  DIVDED 

640  !   BY  N-5.  E.6.  FOR  7  THERMOCOUPLES  Tim  IS  .745,  FOR  8  IT  IS  .85Z 

650  END  IF 

660  NEXT  J 

670  PRINT  "  " 

680  PRINT  "  - 

690  FOR  J-1  TO  N 

700  Tl(J)-A+<B»V1(J))  +  (C»<V1<J)*Z))+(D»<V1(J)-3)>  +  (E»(VKJ)A4)) 

710  TZ(  J)-A+<  B«VZ<  J) )  +  ( C*< VZ( J)*Z)  )  +  (D»<  VZ( J)"3)  )  +  (E«( VZ<  J)  "4)  ) 

7Z0  T3<  J)  =  A+(B»V3<  J)  )  +  <  C»<  V3<  J)*Z)  )+<D»(  V3(  J)*3)  )+<£•<  V3(  J) "4)  ) 

730  T4(  J)=A+<B»V4(  J)  )  +  (  0(  V4<  J  )  *Z  )  )  +  (  D«<  V4<  J  )  *  3  )  >  +  <  E»(  V4(  J  >  "4  )  ) 

740  T5<  J)=A+<B»V5<  J)  )+<  C»<  V5<  J)*Z)  )+<D»<  V5<  J>*  3)  >  +  (E*<  V5<  J)*4)  ) 

750  T6<  J)«A+(B»V6<  J)  >  +  <  C»<  V6<  J)*Z)  >  +  <D«»<  V6<  J)*3>  >  +  <E»<  V6<  J)*4>  ) 

760  T7<  J)«A*(B»V7<  J)  )  +  <  C*(  V7<  J)*Z>  )  +  <D»<  V7<  J)*3>  )  +  <EM  V7<J>"4)  ) 

770  T8(  J)*A+(B«V8<  J  >  )+(  C*<  V8<  J  )  *Z  )  >  +  (  D»<  V8<  J  )  *  3)  )+(  E*<  V8(  J  )  '4)  ) 

780  PRINT  USING  "  1  X ,DDD, 1 X ,4DE , 1 X ,4DE ,  1  X ,4DE . 1 X .4DE . 1 X ,4DE , 1 X ,40E . 1 X ,4DE , 1 X .4 

DE";G< J) ,T1( J) ,TZ( J) ,T3( J)  ,T4( J)  , T5( J),T6(  J ) , T7< J) , T8( J) 

790  NEXT  J 

800  OUTPUT   ©Trail  1 ;G( •  )  ,  T1 ( ♦  )  ,  TZ( ♦ ) , T3( • ) , T4(  • ) , T5< • ) , T6( • ) . T7(  » ) . T8(  • ) 

810  END 
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APPENDIX  B 
EQUATIONS  FOR  DETERMINING  FLUID  PROPERTIES 


P  =  -7.945E-10T^  +  5.736E-5  T FaM  -  9.781E-3 


v  =4.697E-9T*/Uf  -  2.537 E-6  7^  +  1.342E-3 


[i  =  3.235  E-7  T2FEM  -  2147  E-3  TFEM  +  0.03617 


v  =  Vfl 


k  =  lA2E-3TFaM  +0.187 


/V  =  4.66E-3  r*      -  2.92  T__.,  +  463 


Curve  fits  were  determined  for  range  of  temperatures  from  280  to 
300  K. 
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APPENDIX  C 
SAMPLE    CALCULATIONS 

The  calculations  were  based  on  temperatures  derived  from  a  run 
with  all  blocks  powered  to  1.0  watt,  no  shroud  in  place.  When 
appropriate,  values  are  calculated  for  block  number  6. 

A.      CHARACTERISTIC  DIMENSIONS 


Perimeter  = 


=  0.2393  m 


2(0.00787  +  0.02388)  +  (2(2(0.00787  +  0.00610)))' 
.+  (2(2(0.02388  +  0.00610))) 


Al  =  (0.02388X0.00787)  =  L879  x  10~4m2 


4,  =  (0.02388)( 0.00610)  =  1.457  x  10~V 


A2  =  (0.00787X0.00610)  =  4.80  x  10"5m2 


ArorAL=(l.879x  10"4)  +  (2)(l.457  x  10"4)  +(2)U.80x  10"5)  =  5.75  x  10"w2 


Atotaia  =  (5.75x  10"4-L879x  10"4)  =  3.87  x  10~4m2 


r      5.75  x  10         ,  .ftu  irr3 
L=      0.2393      =2-40x10    m 
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B.   SURFACE  AVERAGED  TEMPERATURES 


(19.82)(  1.879  x  10"4)  +  (19.48)(4.80x  10"5) 
+  (19.49)(4.80x  10"5)  +  (20.62)(L457  x  10"4) 
+  (20.26)(L457x  10"4) 


AVG  (8) 


(5.75  x  10"4) 


=  20.08  °C 


( 19.41)(  1.879  x  10"4)  +(18.99)(4.80x  10"5) 
+  (19.05)(4. 80  x  10"5)  +  (19.88)(L457  x  10"4) 
+  (20.08)(L457x  10"4) 


AVG  (7) 


(5.75  x  10"4) 


=  19.64^ 


( 19. 16)(  1.879  x  10"4)  +  (18.67)(4.80x  10"5) 
+  (18.61)(4.80x  10"5)+  (20.00)(L457  x  10"4) 
+  (19.57)(L457x  10"4) 


AVG  (6) 


(5.75  x  10~4) 


=  19.40  °C 


"(18.96)(  1.879  x 

10"4)+(18.33)(4.80x 

lO"5)    1 

+  (18.47)(4.80: 

<  10"5)  +  (19.68)(L457 

x  lO"4) 

_+(19.37)(l457 

x  10-4) 

AVG  (5 ) 


(5.75  x  10"4) 


=  19.16°C 


AVG  (4) 


(18.65)(  1.879  x  10"4)  +(18.30)(4.80x  10"5) 
+  (18.39)(4.80x  10"5)+  (19.75)(L457x  10"4) 

,+  (19.06)(l457x  10"4) . 

(5.75  x  10"4) 


=  18.99  °C 
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(17.76)(4.80x  10"5)  +  (18.10)(4.80x  10"5) 
+  (19.18)(l457  x  10"4)  +  (18.87)(l457  x  10"4) 


AVG  (3) 


(3.87x  10"4) 


=  18.77  °C 


"( 18.07)(  1.879  x 

10"4)  +(17.40)(4.80x 

lO"5)    1 

+  (17.79)U.80> 

<  10_5)  +  (18.85)(L457 

x  lO"4) 

_+(18.56)(l457 

x  lO"4) 

AVG  (2) 


(5.75x  lO-4) 


=  18.32^ 


(17.73)(  1.879  x  10-4)  +  (17.25)(4.  80  x  10-5) 
+  (17.25)(4.80x  10_5)+  (18.60)(l.457  x  10"4) 
+  (17.78)(L457x  10-4) 


AVG  (1) 


(5.75  x  10-4) 


=  17.89  °C 


C.      TFILM 


j  =  20.08+  19.64+  19.40+  19.16  +  18.99+  18.77+  18.32+  17.89  _  {g  Gy>c 


T  19.03+13.23      1fi  nor 

FILM   "  2  " 


D.     TEMPERATURE  EXCESSES 

Using  data  from  block  6,  temperature  excess  is  evaluated  as: 


l\Tf  =  19.16-  13.23=5.93^ 


ATR  =18.67-  13.23=  5.44  <^C 
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ATL  =  18.61-  13.23=  5.38°C 


ATT  =  20.00  -13.23  =6.71°C 


ATa  =  19.57 -13. 23  =  6.34°  C 

D 


ATH  =  24.20  -  13.23  =  10.97  °C 


E.      CONVECTED  HEAT  FLUX 


q„=        0.984  5240h^ 

L879x  10"  m 


F.      TEMPERATURE  SCALING  FACTOR 


(5240)(240xlQ-3)_ 
1     ~  (0.597)  ~ZLl  ^ 


G.     MODIFIED  GRASHOF  NUMBER 


4 

n  „     (9.81)(l63  x  10~6)(  5240)  (  2.40  x  10~3)       ._, 
Gr  *= =  371 

(0.597)(ll2x  10"8) 


H.  HEAT  TRANSFER  COEFFICIENT 


*=         <524°>         =849-" 


(19.40-  13.23)  m2K 
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I.       NUSSELT  NUMBER 


(849)(2.40x  1(T3) 
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